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Dynamical decoupling is a powerful technique for extending the coherence time (T2) of qubits. 
We apply this technique to the electron spin qubit of a single nitrogen-vacancy center in type Ha 
diamond. In a crystal with natural abundance of 13 C nuclear spins, we extend the decoherence time 
up to 2.2 ms. This is close to the Ti value of this NV center (4 ms). Since dynamical decoupling 
must perform well for arbitrary initial conditions, we measured the dependence on the initial state 
and compared the performance of different sequences with respect to initial state dependence and 
robustness to experimental imperfections. 



I. INTRODUCTION 

Sufficiently long coherence times (T 2 ) are amongst the 
most important criteria for the realization of scalable 
quantum information processors. [1] In solid state sys- 
tems, it is often possible to extend the coherence time 
by removing sources of decoherence from the host ma- 
terial, e.g. by isotopic engineering. [2, 3] If this is not 
possible, or to further reduce the detrimental effects of a 
noisy environment, it is possible to use a series of con- 
trol pulses applied to the qubits. This approach, termed 
Dynamical Decoupling (DD), refocuses the interaction 
between system and environment by applying a strobo- 
scopic sequence of inversion pulses to the qubits. It has 
been tested on different systems with environments con- 
sisting of electronic or nuclear spin baths. [4-9] The tech- 
nique allows, e.g., to extend the decoherence time T2 or 
to increase the sensitivity of magnetic field sensors to ac 
magnetic fields. [8, 10] 

The usefulness of this approach for practical applica- 
tions depends on the robustness of the performance with 
respect to experimental imperfections, such as finite pre- 
cision of control field amplitudes. Since the number of 
control operations required for effective DD can be very 
large, it is essential that the errors from the individual 
control operations do not accumulate. This goal can be 
achieved by using fault-tolerant sequences, which are de- 
signed in such a way that the errors of any individual 
control operation is cancelled by all the other operations 
of a cycle. [11-13] 

Here, we test this approach on a single Nitrogen- 
Vacancy color center (NV center) in diamond. The ef- 
fectiveness of the DD technique has been tested with 
different types of diamond; de Lange, et al.[7] used a 
diamond with a relatively high concentration of nitro- 
gen impurities, which generates a strong electron-spin 
bath. The other studies employed CVD grown diamond 
crystals with reduced native nitrogen concentrations, in 
which the nuclear spin bath of the 13 C nuclear spins is 
the major source of decoherence. [8, 11, 14] Here, we ex- 
tend these earlier studies by testing the limits of DD in 
two respects: first, we use DD to extend the dephasing 



time T2 to nearly the value of the spin-relaxation time 
T\. To reach this limit, we have to apply hundreds of re- 
focusing pulses. With such a large number of pulses, the 
unavoidable imperfections of the individual gate opera- 
tions would normally destroy the coherence of the qubits. 
To avoid this, we use robust DD sequences, which were 
designed such that the errors due to the individual gate 
operations do not add up but cancel each other over the 
complete DD cycle. 

When the effectiveness of DD sequences is demon- 
strated, it is often tempting to consider its performance 
for a specific initial condition that is left unperturbed 
by the combination of DD and the environmental per- 
turbation. However, in actual quantum computers, the 
state of the system is generally not known. Accordingly, 
any serious test of DD performance must consider the 
performance of arbitrary initial conditions. We therefore 
present here a study in which we compare the perfor- 
mance of different DD sequences over oil possible input 
states. 



II. EXPERIMENTAL SETUP 

In the experiments described below, we use a home- 
built confocal microscope for optical addressing and de- 
tection of single NV centers. Fig. 1 (a) shows, as an 
example, the scan image of one NV center. 

A diode-pumped solid state laser working at a wave- 
length of 532 nm provides the optical excitation. The 
acousto-optical modulator used for pulse shaping of the 
laser provides an extinction ratio of more than 57 dB. To 
verify that we detect signals from single NV centers, we 
use a standard photon-antibunching configuration with 
single photon detectors in combination with a fast time- 
correlation card. Fig. 1 (b) shows a typical correlation 
signal trace. 

We generate the microwave pulses with an arbitrary 
waveform generator. Working at 2 GS/s, it generates 
the pulses around a carrier frequency of 250 MHz. These 
pulses are upconverted to the target frequency in the 
range of 2.6 - 3.2 GHz by mixing them with a cw signal 
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FIG. 1. (color online)(a) Confocal image (0.8 ^m*0.8 fan) 
of single NV center (b) Second order correlation measure- 
ment by photon antibunching experiment, (c) Ramsey fringe 
(d) Fourier transformation of Ramsey fringe curve shows only 
three peaks due to hyperfine interaction with 14 N in the NV 
center. 

from a stable microwave source. Unwanted sidebands are 
filtered out by properly selected band-pass filters. The 
amplified pulses pass through a 20 /im diameter copper 
wire attached to the surface of the diamond crystal and 
terminated by 50 fl. The leakage of microwave power at 
pulse-off time is < —60 dB. Under typical experimental 
conditions, the duration of a 180 degree pulse is 40 ns. 
A permanent magnet is used to apply a static magnetic 
field, which is roughly aligned with the direction of the 
principal axis of the NV center. A type Ha diamond 
crystal with a N concentration < 5 ppb is the host ma- 
terial of the investigated NV centers. All the measured 
curves in the present work are normalized in order to in- 
dicate the probability (p) of bright state (ms = 0). As a 
reference for the normalization, maximum and minimum 
values from Rabi oscillation curve were taken for p = 1 
and p — respectively. 



III. EXPERIMENTAL RESULTS 

A. Spin system 

In the natural abundance diamond crystal with re- 
duced N concentration, the 13 C nuclear spin form a spin 
bath that induces decoherence on the NV center. [3, 15] 
In the case of Hahn echoes, the deocherence due to the 
13 C nuclear spins is reduced significantly when the free 
precession period between the excitation and refocusing 
pulse is a multiple of the nuclear spin Larmor preces- 
sion period. [16] For the present study, we selected an NV 
center that has no 13 C nuclear spins in its immediate 
vicinity, so the coupling to the nuclear spin bath is rel- 



atively weak and the decoherence correspondingly slow. 
Figure 1 (c) shows the Ramsey fringe signal measured 
from this center and Figure 1 (d) its Fourier transform. 
We applied a magnetic field of 6.8 mT along the direction 
of the symmetry axis and excited the transition at 2.68 
GHz with resonant microwave pulses. 



B. Dynamical decoupling 

The oldest and simplest pulse sequence for dynamical 
decoupling is the CPMG sequence - a train of equidistant 
7r-pulses around the same axis. In Fig. 2, we compare 
the signal decays for different CPMG sequences, starting 
with the Hahn echo at the bottom and experiments with 
increasing numbers of refocusing pulses towards the top. 
In each curve, we observe an initial fast decay, which is 
followed by a series of revivals. The revivals are sepa- 
rated by 2Ntl, where N is the number of it pulses in the 
sequence and tl ~ 73 ^s is the Larmor period of the 13 C 
nuclear spins. This implies that the refocusing is most 
effective when the nuclear spins constituting the envi- 
ronment have completed an integral number of Larmor 
precessions between the pulses. The dashed curves in 
Fig. 2 represent functions exp(— (y-) n ) with decay times 
Ti and exponents n fitted to the maxima of the Larmor 
revivals in the experimental data. The decoherence times 
Ti are also noted in the figure legend. 

Figure 3 summarizes the increase of the decoherence 
time with the number of refocusing pulses. The deco- 
herence times from the CPMG experiment are shown as 
black squares, with filled squares representing the values 
obtained by fitting the decay of the Larmor revivals (the 
dashed curves in Fig. 2), while the open squares represent 
the fit to the initial decay. In both cases, the decay time 
increases with the number of pulses in the DD sequence. 
If we take the values from the maxima of the Larmor 
revivals, they increase from rs 0.3 ms for the Hahn echo 
up to 2.2 (±0.2) ms for CPMG 128, which is the longest 
value ever reported for the NV center. 

Apart from the increase of the decoherence time with 
the number of pulses, we find that also the shape of the 
decay curve changes. In Fig. 3, we quantify this by the 
exponent n, which changes from ss 3 for the Hahn echo to 
« 1 for CPMG 128. This result indicates that different 
processes are responsible for the observed decoherence. 
For relatively long pulse separations, as in the case of the 
Hahn echo, the dominant environmental interaction is 
the hyperfine coupling to the nuclear spins. The shortest 
relevant timescale for the nuclear spins is the Larmor 
precession, which has a period of tl sa 73 fis. Since this 
evolution is coherent, the associated decay is refocused 
at specific times, giving rise to the Larmor revivals [17]. 
On a slower timescale, the nuclear spins undergo mutual 
spin flips, which also contribute to the decoherence of the 
electron spin. 

When the spacing between the pulses becomes short 
compared to the Larmor period, the refocusing of the de- 
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FIG. 2. (color online) Coherence decay curves for different 
dynamical decoupling sequences. The bottom curve shows 
the result from the Hahn echo, the higher curves correspond 
to the CPMG sequence with increasing number of refocusing 
pulses. The gray dashed lines represent fits to the function 
exp( — (t^-)"). The initial state was chosen along the rota- 
tion axis of the refocusing pulses. The top curve shows Ti 
relaxation of bright state. (Ti = 4 ms) 



coherence due to the nuclear spin bath becomes almost 
perfect. This is demonstrated in Fig. 3 by the conver- 
gence of the curves from the decay time of the envelope 
with that of the initial decay. Both time scales approach 
a limiting value of ss 2.2 ms, which is only slightly shorter 
than the longitudinal relaxation time for this spin, which 
is Ti ~ 4 ms. Apparently, the remaining decoherence 
processes causing this decay are no longer due to the 
interaction with the nuclear spins, but arise from pro- 
cesses with significantly shorter correlation times, which 
cannot be refocused by dynamical decoupling. Possible 
candidates for these processes include interactions with 
paramagnetic centers, electric fields [18] or phononic pro- 
cesses. This interpretation is compatible with the obser- 
vation of a purely exponential decay, which is expected 
if the correlation time of environment is shorter than in- 
terpulse delay of dynamical decoupling. 



C. Parallel and orthogonal initial states 

The above experiments were performed with one spe- 
cific initial condition, for which the CPMG sequence is 
optimized. However, in quantum information processing, 
the initial state is in general unknown, and it is impor- 
tant to use a pulse sequence that performs well for ar- 
bitrary initial conditions. To evaluate the effectiveness 
for arbitrary initial conditions, we now consider the ef- 
fect of the DD sequence on two orthogonal initial con- 
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FIG. 3. (color online) T2 and T2 as a function of the number 
of pulses in CPMG and KDD sequences. Black symbols are 
CPMG and red KDD. Solid symbols represent T2S obtained 
from the envelope of the revivals, and open symbols are T°s 
from the initial decay. Green symbols represent the exponent 
n in the decay function exp(— (-jr) n ). The initial state was 
chosen along the rotation axis of the CPMG pulses. 



ditions \x) = (|0) + and \y) = (|0) +i|l))/ v / 2. 

Fig. 4 shows the decay of these states for the DD se- 
quences CPMG, KDD, and XY-4 as a function of the 
number of refocusing pulses. In the case of the CPMG 
sequence, we observe a strong asymmetry between the 
two initial condition: If the initial spin polarization is 
oriented along the rotation axis of the refocusing pulses, 
the state preservation works extremely well. If the ini- 
tial state is perpendicular to the direction of the pulses, 
the first w 10 pulses destroy the coherence almost com- 
pletely. This effect is expected if the rotation angle of 
the pulses deviates from the ideal value of 7r.[9] Such a 
flip angle error has no effect on a state that is polarized 
along the rotation axis, but for orthogonal orientations, 
the flip angle errors of the individual pulses accumulate 
and cause an unwanted evolution of the system. 

Such a strong dependence on the initial state of the 
system is not compatible with quantum information pro- 
cessing, where the state is not known in general. Univer- 
sal DD sequences must perform well for arbitrary quan- 
tum states and they have to be robust against pulse er- 
rors. The earliest sequence that fulfills this requirement 
is the XY-4 sequence [19, 20] and its derivatives that in- 
crease the robustness by combining different versions of 
the basic XY-4 cycle into longer cycles with better com- 
pensation of pulse imperfections [21, 22]. A more recent 
sequence that shows even better performance is the KDD 
sequence [13], which is based on an expansion scheme 
developed by Cho et al [23]. As shown in Fig. 4, the 
performance of both sequences does not depend on the 
initial state, within the experimental uncertainties. The 
performance of the XY-4 sequence is comparable to that 
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FIG. 4. (color online) Survival probabilities of two orthogo- 
nal initial states, |sc) (solid symbols) and \y) (empty) during 
different sequences of equidistant refocusing pulses for CPMG 
(red), KDD (blue), and XY-4 (green). The time interval be- 
tween the refocusing pulses was fixed to 0.8 ps. 
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of CPMG_y for < 20 pulses but deteriorates for longer 
sequences, while the KDD sequence matches the perfor- 
mance of CPMG_y for both initial conditions, indicating 
that this sequence largely compensates the flip angle er- 
rors of the individual pulses. 



FIG. 5. (color online) Experimental survival probability for 
the sequences CPMG, KDD and XY-4 after 20 DD pulses 
as a function of the parameters 9 and <j> that define the ini- 
tial state of eq. (I). Top left corner illustrates Bloch sphere 
representation of arbitrary state of NV center. 



E. Effect of pulse imperfections 



D. Arbitrary initial states 

Fig. 5 extends this study to arbitrary initial conditions: 
for these experiments, we prepared the initial states 

|^)=cos(|)|0)+sin(|)e^|l), (1) 

where the angles and <fi, which parametrize the state, 
correspond to spherical coordinates on the Bloch sphere. 
Fig. 5 shows the observed survival probability |(^|^j)| 2 , 
with \tpp) representing the state of the qubit after 20 
refocusing pulses. The three panels represent the sur- 
vival probability as a function of 6 and cf> for the three 
DD sequences CPMG, XY-4 and KDD. All three se- 
quences reach almost perfect survival probabilities for 
9 = <j) = 7r/2, which corresponds to the y initial condition 
where the spins are aligned with the rotation axis of the 
CPMG pulses. In the case of CPMG, the performance 
deteriorates when the initial condition deviates signifi- 
cantly from this situation. This behavior is an indication 
that the decay of the system fidelity is dominated not 
by the environment, but rather by pulse imperfections 
[9,13]. 



Apart from finite pulse durations, the two dominant 
experimental imperfections that limit the performance of 
experimental dynamical decoupling are flip angle errors 
and offset errors. In Fig. 6, we experimentally investi- 
gate the effect of these imperfections on the performance 
of different dynamical decoupling sequences. In the left 
hand columns, we measured the remaining spin polariza- 
tion after 20 and 200 pulses as a function of the actual 
flip angle of each pulse. For 20 refocusing pulses and lon- 
gitudinal initial conditions, flip angle errors up to ± 10 
% do not reduce the fidelity. For higher errors, we find a 
bi-quadratic reduction (not shown here). This is consis- 
tent with theoretical considerations. [24] For the robust 
sequences, whose performance does not depend signifi- 
cantly on the initial condition, KDD clearly outperforms 
XY-4. This is particularly well visible for 200 pulses. In- 
terestingly, XY-4 shows strong oscillations as a function 
of the flip angle error. This implies that XY-4 becomes 
very sensitive to flip angle errors if more than 100 refocus- 
ing pulses are used, which limits its usability. The right 
hand column shows the measured behavior as a function 
of the offset error for 20 refocusing pulses. For all three 
sequences, we observe significant reductions of the fidelity 
if the offset exceeds ps 2 MHz. This is significantly less 
than the Rabi frequency (12.5 MHz) of the pulses used 
for these experiments and comparable to inherent offsets 
like the hyperfine interaction with the nitrogen nuclear 
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FIG. 6. (color online) Left column : Survival probability of 
the electron spin state as a function of the fractional flip angle 
error. The initial state was chosen along the rotation axis of 
the CPMG pulses. The numbers in the parentheses represent 
number of refocusing pulses in DD sequence. Right column 
: variation of survival probability as a function of frequency 
offset errors on pulses. The numbers in the parentheses rep- 
resent number of refocusing pulses in the DD sequence. 



spin. 

In order to better understand the processes that re- 
duce the process fidelity in these cases, we performed 
state tomography after 20 refocusing pulses for the ini- 
tial state along the x-axis (0 = 5, <j> = 0). From the state 
tomography, we could estimate a flip angle error of 1.6 
%, which corresponds to 2.9 degree deviation from the 
targeted 180 degree rotation. This value is quite consis- 
tent with our hardware limit on the width of time step. 
The arbitrary wave form generator used for pulse gener- 
ation has 0.5 ns step width (2 GS/s). Given the duration 
of the refocusing pulse of 40 ns, the digital time resolu- 
tion of the pulse results in a flip angle error of 1.2 %, 
which is not far from value given above and does not 
yet include effects due to amplifier nonlinearities or tran- 
sient effects. Frequency offset error was also observed, 
and this was - 4 % (0.5 MHz) of the Rabi oscillation fre- 
quency (12.5 MHz). This amount of estimation can be 



supported by the results of the frequency offset errors in 
Fig 6. The central positions of all the three curves are 
slightly shifted to positive frequency, and this deviation 
is around 0.5 MHz. The origin of this frequency offset 
error is unclear to us, because we tune the frequency of 
the microwave pulses to the center of the spectrum mea- 
sured by the Ramsey fringe experiment in Fig. 1 (d) as 
precisely as possible, expecting that the off-resonance ef- 
fect of the two outer hyperfine lines cancel each other. 
There is one way to confirm whether the off-resonance 
effect generates frequency offset errors on the pulses: do- 
ing experiments at the excited state level anti-crossing 
point (514 G), only one of the 14 N nuclear spin states 
will be populated, [25] and this will make it available to 
remove off-resonance effect completely. 



IV. CONCLUSIONS 

Dynamical decoupling is an effective method for in- 
creasing the coherence time of quantum bits, such as spin 
qubits in solids. As long as the fluctuation of the noise 
sources is comparable to or slower than the system's co- 
herence time, DD can increase the survival time of the 
quantum information. With CPMG sequences, we have 
shown that T2 of a single NV center in natural abundance 
diamond can be extended up to 2.2 ms by strongly sup- 
pressing the fluctuating fields from the 13 C nuclear spin 
bath. This is close to the Ti limit, which is determined by 
different, more rapidly fluctuating processes. We found 
that the recently introduced DD sequence KDD is re- 
markably robust and works for arbitrary initial condi- 
tions. It is much less susceptible to experimental uncer- 
tainties, such as flip-angle errors and frequency offsets, 
than the simpler sequences CPMG and XY-4. We ex- 
pect that this robust DD sequence will be used for many 
other purposes requiring long coherence time, such as ac- 
magnetometry or to increase the lifetime of multipartite 
systems, including entangled states. 
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